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Abstract 

We are interested in a nonlocal conservation law which describes the morphodynamics of sand 
dunes sheared by a fluid flow, recently proposed by Andrew C. Fowler and studied by d [21 - We 
prove that constant solutions of Fowler's equation are non-linearly unstable. We also illustrate this 
fact using a finite difference scheme. 
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1 Introduction 

Partial Differential Equations with nonlocal or fractional operators are widely used to model scientific 
problems in mechanics, physics, signal processing and other subjects, see for example [3 ] and references 
therein. We consider in this paper a nonlocal conservation law which appears in the formation and 
dynamics of sand structures such as dunes and ripples 0QTI. Namely, Fowler (GHHm) introduced the 
following equation 

(d t u(t,x) + d x ^(t,x)+X[u(t,-)]{x)-d^ x u(t,x) = t€(0,T),xeR, 
\u(0, x) = uq{x) x € R, 

where u = u(t, x) represents the dune height, uq G L 2 (R) is an initial condition, T is any given positive 
time and I is a nonlocal operator defined as follows: for any Schwartz function <p G <S(R) and any 

T[<p]{x) := / |CrV(z-CR- (2) 
J o 

Equation CD) is valid for a river flow over an erodible bottom u(t, x) with slow variation. The nonlocal 
term appears after a subtle modeling of the basal shear stress. This operator appears also in the work of 
Kouakou & Lagree ifTOlfTTI 

The nonlocal term I can be seen as a fractional power of order 2/3 of the Laplacian with the bad sign. 
Indeed, it has been proved HI 

t (xm - ip") (0 = MO^viO 0) 
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where 



Mo=^ 2 e - ax \^ +i b x m s, 



(4) 



with ax, bj positive constants and T denotes the Fourier transform. One simple way to establish this fact 
is the derivation of a new formula for the operator X 12 

2M(*) = C X f dz, (5) 

with C x = |. 

The operator Z[u] is a weighted mean of second derivatives of u with the bad sign and has therefore an 
anti-diffusive effect and creates instabilities which are controlled by the diffusive operator —d 2 x . This 
remarkable feature enabled to apply this model for signal processing. Indeed, the diffusion is used to 
reduce the noise whereas the nonlocal anti-diffusion is used to enhance the contrast Q. 



Remark 1. For causal functions ( i.e. <p(x) = Ofor x < 0), this operator is up to multiplicative constant, 
the Riemann-Liouville fractional derivative operator which is defined as follows / I72l/ 



i r*°v>"{x-$ d-^ d^ 



r(2/3)7 |e|V3 s dx-W w dx 4 / 3 

where V denotes the Euler function. 

Recently, some results regarding this equation have been obtained, namely, existence of travelling- 
waves u^(t,x) = (j){x — ct) where <p G C^(R) and c G R represents wave velocity, the global well- 
posedness for L? -initial data, the failure of the maximum principle, the local-in-time well-posedness in 
a subspace of C\ UlEl and the global well-posedness in a L 2 -neighbourhood of C\, solving namely 
for u = Uff, + v, with v G L 2 (R) 0. For this purpose, the following Cauchy problem solved by the 
perturbation v has been considered: 

fd t v(t, x) + d x (^ + U(j) v)(t, x) + X[v(t, -)}{x) - d 2 xx v(t, x) = t G (0, T), x G R, 
vq(x) x G R, 



\v(0,x) 



where vq G L 2 (R) is an initial perturbation and T is any given positive time. 

Let us note also that any constant is solution of the Fowler equation. We shall prove below that these 
solutions are unstable. 

To prove that a solution is unstable when 4> is constant, we introduce the notion of mild solution (see 
Definition [B based on Duhamel's formula ©. We also give some numerical results that illustrate this 
fact. 

The remaining of this paper is organized as follows: in the next section, we define the notion of mild 
solution and we give some results. Section |3]contains the proof of the instability. We introduce in section 
H]an explicit finite difference scheme for which we give numerical simulations to illustrate the theory of 
the previous section. 
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2 Duhamel formula and some results 



Definition 1. Let be a constant, T > and v G L 2 {R). We say that v G L°°((0, T); L 2 (R)) is a 
mild solution to (Q if for any t G (0, T), 



x) = K(t, •) * uo(a:) - / d x K{t - s , •)*( — ) ( s , .)( x ) (8) 



o 



where K(t, x) = T 1 (e ^v)) (x) is f/ie kernel of the operator X — d 2 x + U0<9 X and 4>x(0 = 47r 2 £ 2 — 

The expression ([8]) is the Duhamel formula and is obtained using the spatial Fourier transform. The 
use of this formula allows to prove the local-in-time existence with the help of the contracting fixed point 
theorem. The global existence is obtained thanks to L 2 a priori estimate. We refer to (Hill for the proof. 

Lemma 1. Let t > 0. Then, d x K(t, •) G L 2 (R) and satisfies 

\\d x K{t,-)\\ L 2 {m <c(t-V A + e at ), (9) 

where C is a positive constant independent oft. 
Proof. By Plancherel formula, we have 

l^(t,-)ii! 2(R) = ii^(^(t,-))ii! 2(K ) = iie^2^ e -^)|ii 2(s 

r+oo 
Jo 



Let £ > such that for all £ > £ > 

4ir 2 e-ax£ i/3 >e. 
Then, if we denote by a = — mini2e(0x) (see Figure©, we have 

f£o , /"+00 



\d x K(t,-)\\ 2 L2{R) < Svr 2 y e 2 e- 2 *( 4 ^ 2 -^ 4/3 )d£ + 8vr 2 y ^-a*^ 8 -^ 8 ) ^ 



< Svr 2 / ? °e 2 e 2to de + 87r 2 /^"Ve" 2 * 3 d£, 

Jo Jo 

t3 /|_2 /*-)-oo 

< 87r 2k e 2to + l^i-3/2 / e 2 e -^ de> 

3 y2 Jo 

< C -f t -3/2 + e 2aA j 



which completes the proof of this lemma. ■ 

Remark 2. Using again Plancherel formula, we have for any initial data vq with values in L 2 (JH), the 
following L 2 -estimate £/]/.• 

IKV)||l2(k) < e at \\v \\ L 2 (R) . 
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3 Instability 

In this section, we give the proof of the instability. 

Let us first note that a fundamental property of the kernel K is the non-positivity [T), see Figure Q] 
This feature enabled to prove the failure of the maximum principle. We use again this property to show 
that the constant solutions of the Fowler equation are unstable. 

1.2 
1 
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Figure 1: Evolution of the kernel K for t = 0.1 and t = 0.5 s 



Theorem 1. Any constant solution of the Fowler equation is unstable, i.e. 3e > 0, V5 > 0, 3vo 6 
L 2 (M) an initial data with | |f o| lz, 2 (M) — <5 an d *o > such that \\v(to, -)||l 2 (r) > £ > where v(t, ■) is 
the solution of © with v(0, ■) = vo- 

Proof. We adapt the strategy used by De Bouard in making it more elementary in our setting. 

We denote by S(t) the linear semigroup associated with equation © i.e. S(t)w = K(t, •) * w. Hence, 

we have 

1 /"* 

T(t)v := v(t, •) = S(t)v - - / d x K{t - s, •) * w 2 (s, •) ds, (10) 



2 „ 

where T denotes the nonlinear semigroup associated with equation ©. 
First, by Young inequality, Lemma[T]and Remark[2l we have 



1 



2 

< C - 
~ 2 



\\T(t)v - S(t)v \\ L 2 {R) < - \\d x K(t- s,-)\\ L2{R) \\v(s,-)\\ L 2 m ds 



o 



o 



( t _ a )-3/4 + e «(t-.) e tt"|| U0 ||* ds 



,2asii„, ||2 



< Ht)\\vo\\h {R) , (11) 

where = § e 2 "*^ 1 / 4 + ±e a *) > for all t G (0, T). 
Hence, we obtain for to > fixed, 

||T(i )uo - <S , (to)uo|| < & olbo|li2 (M) , (12) 

where bo = b(to). 

Assume now that the constant is a stable solution of the Fowler equation i.e, 

Ve > 0, 3/7 > 0,||«o|| < V IN*, -)IIl 2 (r) < e > 
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for all positive t. Let then 

ato _ i 

0<£< 



16&o 

and let r\ > be such that ||uo||z,2( K ) < rj =4> \\v(t, -)|| < e, for all positive t. 
Let N be an integer large enough such that 

e at N > f 1_ 

- 4r/6 

Let us now consider the following initial data: vq = <5^o> where J-"(tt;o) = ^ e l[c,rf]' with < c < d 
satisfying (see Figure [2]): 

• Re (4>x) (c) < Re (<fc) (d) = -/3 < 0, 
. for all £e M,Re(<fc)(£) < -ft 

• j3 = a — 7 where < 7 < 
And 5 is defined as 

at _ 1 

^ = e- QjVf ° - . (14) 
Ab 



The shape of wo is given by 

wo (a;) 



> sin( ^ x) -^ if x ^0, 



y/d — c I d — c otherwise. 
It follows from (fT3l) and (fT4l that ||vo||l 2 (R) = ^ ^ >7> an d thus for this initial data, we get that 




Figure 2: Behaviour of Re (</>j) 

IKV)IIl2(r) < e for all t > 0. 
For n G {0, ■ ■ ■ , iV}, we have 

n-l 

u(nt ) = S{nt )v + Y,S{(n-l-k)t )[v((k + l)t )-S(t )v(kt )}, 

n-l 

= S{nt )v + Y,S((n-l- k)to) [T(t )v(kt ) - S(t )v(kt )} . 

k=0 
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Hence, we get using (fl2l) and the L 2 -estimate given in Remark[2] 

n-l 

\\v(nt ) - S(nt )v Q \\ LH K) < E e ato{n - k ~ 1) M(k + l)t ) - S(t )v(kt )\\ L 2 { 

k=0 

n-l 



k=0 
n-l 



fc=0 



„anto -i 

^ ato(n-l)e i 

e ai _ 1 



" 4 ' 



Moreover, by Plancherel formula, we have 
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\\ S ( t ) v o\\mR) = W K ( t >-)*Vo\\ L 2 {R) = \\T(K(t,-) *V )\\ L 2(jgr ) , 

\T(K(t,-))(OHvom\ 2 dt, 



d j!_ e -2tRe(fe)(0 £ ^ ) 
a — c 



We finally infer that 



lk(^ r *o)||i2(R) > \\S{Nt Q )v Q \\ L 2^)-\\v{Nt Q )- S{Nt )v Q \\ L 2(^), 

> Se PNU> _ IjxtoN 
4 

e a*0 _ J 

166 

which gives us a contradiction and completes the proof of this theorem. 

■ 

Remark 3. We can give a physical interpretation of this result: aflat profile us = constant is unstable 
under the morphodynamics described by the Fowler model. 

4 Numerical simulations 

The spatial discretization is given by a set of points Xj\ j = 1,...,N and the discretization in time is 
represented by a sequence of times t° = < ... < t n < ... < T. For the sake of simplicity we will 
assume constant step size 5x and 5t in space and time, respectively. The discrete solution at a point will 
be represented by u r - ps u(t n , Xj). The schemes consist in computing approximate values of solution 
to CD on [nSt, (n + l)5t[x [j5x, (j + l)5x[ for n € N and j G N thanks to the following relation: 

A + 2 fa fa* ]j " °' (15) 
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where I$ x is the discretization of the nonlocal term X. We use a basic quadrature rule to approximate X 
given by © 

+00 

XsMj = fe~ 4 / 3 J^r 1 / 3 (<pj-i+i - %Pi-x + ipj-i-x) ■ (16) 
1=1 

Let us note that for compactly supported initial datum the sum is a finite sum. In this section, we want 
to simulate the instability stated previously. To this aim we must be careful to distinguish between the 
instability which stems from the model and numerical instabilities, which could be caused by a careless 
discretization. 

The numerical stability ensures that the difference between the approximate solution and the exact solu- 
tion remains bounded for increasing t for Sx, St given. 

As noticed, the Fowler model amplifies slowly the low frequencies whereas the high frequencies are 
quickly dampened. Thus, the classical notion of ^4-stability (strong stability) and C-stability are not suit- 
able nor desirable. A new definition of stability for this model has been considered in and numerical 
stability criteria have been obtained. 

In the following numerical test, we have to take care to choose Sx and St accurately following the 
stability condition, see [51). We expose in Figure [3(a)] the evolution of an initial flat bottom disturbed by a 
small bump for different times. As we can see the perturbation creates small ripples that grow with time 
without bound. This result proves that the perturbed solution goes away from the non-perturbed solution. 
This illustrates the instability of constant solutions for the Fowler model £[]). Figure [3(b)] shows the log- 
plot for the L 2 -norm evolution of this perturbed bottom, which confirms the validity of the estimate given 
in Remark |2] 




* t(ms) 
(a) Evolution of the solution disturbed. (b) Evolution of log(||i)(t, •) I II 2 ) 



Figure 3: Evolution of a flat bottom disturbed by a small bump. 
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